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ROBUSTNESS OF FLEXIBLE SYSTEMS 
WITH COMPONENT-LEVEL 
UNCERTAINTIES 

Peiman G. Maghami * 

NASA Langley Research Center, Hampton , VA 23681 

Robustness of flexible systems in the presence of model uncertainties at the 
component level is considered. Specifically, an approach for formulating robustness 
of flexible systems in the presence of frequency and damping uncertainties at the 
component level is presented. The synthesis of the components is based on a 
modifications of a controls-based algorithm for component mode synthesis. The 
formulation deals first with robustness of synthesized flexible systems. It is then 
extended to deal with global (non-synthesized ) dynamic models with component-level 
uncertainties by projecting uncertainties from component levels to system level. A 
numerical example involving a two-dimensional simulated docking problem is worked 
out to demonstrate the feasibility of the proposed approach. 


Introduction 

Component mode synthesis is a well established 
area in modeling of flexible systems 1 There are 
several advantages of modeling via component mode 
synthesis. These include: (1) It allows much inde- 
pendence in the design and analysis of components , 
i.e., the work can be performed by different groups or 
organizations. A good example is the International 
Space Station which are built by different companies 
in several countries. (2) It increases the power 
of existing finite-element analysis and design, by 
allowing it to handle excessively large problems; (3) 
It allows a direct synthesis of substructure test data. 
This is particularly useful for very large structural 
systems (such as the International Space Station) 
that cannot be tested as a whole. It can also be 
used as a part of an experimental verification tool 
for substructures before deployment as a connected 
structure, for example, to provide uncertainty mod- 
els for the component dynamics. 

In applications of robust control to flexible sys- 
tems, uncertainty models have been typically as- 
sumed or defined for the combined (full-up) system. 
This approach is acceptable so long as these un- 
certainties (in terms of structures and levels) have 
been verified through experimental work. However, 
this is not the case in many situations. A good 
example is where ground testing of the assembled 


* Senior Research Engineer, Guidance and Control Branch, 
Senior Member AIAA. 

Copyright © 2000 by the American Institute of Aeronautics and 

Astronautics, Inc. No copyright is asserted in the United States under Title 
17 U.S. Code. The U S. Government has a royalty-free license to exercise 
all’ rights under the copyright claimed herein for Governmental Purposes 
All other rights are reserved by the copyright owner. 


structure is not possible, i.e., only component-based 
testing is feasible. In such cases, reliable uncertainty 
models can only be defined for the components. 
Consequently, methods are needed to provide for 
robust design and analysis of flexible systems with 
the uncertainty defined at the component levels. 

In this paper, robustness of flexible systems with 
component level uncertainties are addressed. Fiist, 
robustness of a synthesized system to uncertainties 
in the frequency and damping of its components 
is considered. The approach used for component 
synthesis is based on a slight modification of the 
control-based synthesis algorithm by Maghami and 
Lim. 7 The synthesized plant in transformed coordi- 
nates is placed in a Linear Fractional Transformation 
(LFT) form for robust design and analysis. Next, 
a formulation is also presented through which the 
uncertainties at the component level can be ex- 
tended to system level uncertainties for design or 
analysis. This is quite useful in situations wherein 
while a global model of the assembled system is avail- 
able, uncertainties are defined for the components 
only. The proposed approach is applied to a two- 
dimensional problem involving the docking of two 
components in space. 

Substructure Synthesis 

Assume that two flexible systems are connected 
as shown in Fig. 1. These systems may represent 
two components of a larger flexible body (such as 
spacecraft or aircraft), or they may represent two 
self-sufficient systems that have joined for a specific 
purpose, such as two spacecraft docking. Following 
the approach of Maghami and Lim, 7 the interface 
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Component Breakdown 


Fig. 1 Diagram of the Two-Component System 



Fig. 2 Block Diagram of the Two-Component 
System 

forces and moments (between the two components) 
may be considered as control input forces. A block 
diagram for the system is illustrated in Fig. 2. 
Here, each of the component plants are represented 
by Pi and P 2 , respectively. Associated with each 
component model, there is uncertainty in the Linear 
Fractional Transformation (LFT) form, representing 
parametric and nonparametric uncertainties in the 
component models. The controller, K s , denotes the 
controller that synthesizes the combined dynamics 
of two components in a joined configuration. The 
state-space representation of the dynamics of the 
components is 


component 1 


±\ = A\X\ + Biu + Hipi + E,wi 

(1) 

*/i - C\X\ 
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component 2 

X o = A 2 X 2 + B 2 U + H 2 P 2 + E'l w 2 (5) 

JJ2 = C 2 X 2 ( 6 ) 

Z 2 = F 2 X - 2 (?) 

Vp2 “ L 2 x 2 + D u 2 u + D P2 P2 (8) 


where A t , i = 1, 2 denote the state matrices for each 
of the components. Matrices Bi, Hi, E { , Ci, F { , 
and Li, i ~ 1,2, represent the influence matrices for 
interface inputs, non- interface inputs, uncertainty 
outputs, interface outputs, uncertainty inputs, and 
non-interface output for each component, respective- 
ly. The matrices D Ul , D Pl , D U2 , and D P2 , represent 
feedthrough matrices associated with non-interface 
outputs of each component; xi and x 2 denote the s- 
tate vectors of the components; y\ and y 2 denote the 
interface output vectors, including the displacements 
and velocities at the interface from each component; 
y pi and y P2 denote the non-interface output vectors 
components which includes performance outputs as 
well as measurement outputs for local feedback 
control; z\ and z 2 denote the uncertainty input 
vectors; pi and p 2 represent the non-interface inputs 
which include exogenous disturbances as well as 
control inputs for local feedback control; and w\ 
and W 2 denote the uncertainty output vector. The 
influence matrices E t and Fi,i = 1,2 corresponding 
to component uncertainties will be defined later. 

Now, the derivation in Ref/ is followed to obtain 
the synthesized dynamics of the combined system. 
However, in this paper, the interface output vector 
includes the interface velocity vector as well as the 
displacement vector at the interface. The motivation 
behind including the displacement vectors is to avoid 
the superfluous rigid- body mode that may come 
about from the previous formulation. Let N c denote 
an orthonormal basis for the right null space of 
[ Ci -C 2 ], i.e., 

[ Ci -C 2 } N c = 0 (9) 

and let R c denote an orthonormal complement to 
N c - Transform the combined system via an orthog- 
onal similarity transformation, such that 

(,o) 

Then, the dynamics of the combined system in 
transformed coordinates becomes 

a = A\a + H\p -f E\tu (11) 
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Model Uncertainty 

The synthesized dynamics of the combined system 
is illustrated in Fig. 3. Here, the plant P, 
representing the system dynamics, is given as 


(21) 


The uncertainty block, A include uncertainties from 
both components, i.e., 

Aj 0 
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(22) 


The uncertainties represented by A may include 
parametric and nonparametric uncertainties in each 
of the components. Parametric uncertainties in- 
clude uncertainties in modal frequencies, damping 
ratios, and mode shapes of the components. The 
nonparametric uncertainties considered here are the 
unmodeled dynamics of the components, which 
would typically be due to the truncated modes of the 
components. In this paper, parametric uncertainties 
in the component models in the form of frequency 
and damping uncertainties are considered. 



Fig. 3 Block Diagram of the Synthesized System 

Uncertainty in Frequency and Damping 

If, as typically done, a modal model of the 
components are used to describe their dynamics, 
then uncertainties in modal frequency and damping 
are embedded in the elements of the component state 
matrices .4 1 and .4 2 ■ Consider the first component, 
and assume that there are ny distinct uncertainties 
associated with the Ay. Then, one may write 

.4, Ay +J2 Qi*i (23) 

i— 1 

where matrix Q t represents the influence matrix for 
the Si uncertainty. Now, the influence matrices E { 
and F\ (see Eqs. (1) and (3)) may be obtained from 
decomposition of matrices Qi,i ~ 1,... ,n i- Factor 
each matrix as 

Q z = E u Fu (24) 

where Eh £R n E u Gft r ‘ x ", and r t denotes the 
rank of Qi> Singular Value Decomposition (SVD) is 
a good way of performing this factorization. Then, 
the influence matrices E\ and F\ may be written as 

E\ = [ E \ i E\2 ••• E \ n i ] (25) 


F x = 


Fxx 

F\2 


Fin* J 


(26) 


With matrices Ex and F\ defined, the synthesized 
system represented by Eq. (21) is in suitable form 
for robustness analysis or control design. 

Component Uncertainty at System Level 

In some applications, it may be the case that 
an overall model of the combined structure is 
available either from full-up structural analysis or 
system identification, i.e., no component synthesis 
is required to generate the model. However, the 
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model uncertainty is defined at the component level, 
where substantially more effort is typically spent 
for modeling, testing, and verification. In order to 
perform robustness analysis or robust control design, 
the uncertainty formulation has to redefined at the 
system level. Assume, without any loss of generality, 
that the overall model at the system level is in a real 
modal form. Then, the model of the system may be 
written as 



A s H s E s 
N s D s 0 
F s 0 0 


Here, A s has the following form 

0 . . . 0 




0 


with 


As, 



0 

Ain 

-An 

A„ 

A ift 


x s 

P 


(27) 


(28) 


(29) 


where X lR and A 2/ denote the real and imaginary 
part of the ith eigenvalue pair of the system. Here, 
we assume that the eigenvalues are ordered from 
smallest to the largest. From component synthesis, 
we know r that the first m eigenvalues of matrix A i 
should closely match the eigenvalues of the A s . m 
is typically chosen as the number of system modes 
that are within half the bandwidth of the component 
dynamics. Let matrix \ denote a real similarity 
transformation through w T hich matrix A\ is put in 
real modal form, with eigenvalue blocks in increasing 
order, and let Such a transformation may 

be obtained from the complex eigenvectors of A\ by 
collecting the real and imaginary parts of the eigen- 
vectors corresponding to complex conjugate pairs of 
eigenvalues. Applying this similarity transformation 
to the system in Eq. (21) gives 
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(30) 

Comparing the system in Eq. (30) with that given 
by Eq. (27), and keeping in mind that the first m 
eigenvalues of the tw r o systems should match, it is 
observed that the uncertainties at the component 
level can not be put exactly in the same structure 
at the system level, mainly due to the additional 
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15 m 5 m 

Component No. 1 Component No. 2 

Fig. 4 Components of the Free-Free Problem 


dynamics in the synthesized model that is neglected. 
However, a reasonable approximation may be made 
by choosing E s and F s as 


E s - $iEi ; F s — F\x\ (31) 


where 4>i and denote partitions formed from 
the first 2 x m rows and 2 x m columns of $ 
and x, respectively. With E s arid F s defined, the 
global system is prepared for robust design and 
analvsis. It should be noted that if more than 
one rigid-body mode are present in the system 
(e.g., space applications), then it w r ould be unlikely 
that the a rigid-body mode in A\ would have a 
exact corresponding mode in A s . This is because 
the eigenvectors corresponding to rigid-body modes 
(with repeated eigenvalues at zero) are not unique. 
One way to handle this is to replace the partition 
in matrices A s , H s , and N s , corresponding to the 
rigid-body modes, with those of \Pj4ix, ^ H ] , and 
NiX, respectively. This note also applies to repeated 
flexible modes, except in many of such cases unique 
eigenvectors may exist. 


Numerical Results 

The proposed approach for robust design and 
analysis under component-level uncertainty models 
is applied to a two-component planar problem shown 
in Figure 4. Here, a docking scenario involving two 
flexible space systems is considered. The tw r o sys- 
tems are modeled as planar free-free Euler-Bernoulli 
beams. The material and geometric properties of 
the system were chosen to provide considerable 
modal content in the low T -mid frequency range, to 
make the synthesis and control design task more 
challenging, and are provided in Table 1. It is 
assumed that in the docked configuration component 
no. 1 would provide attitude and position control for 
the combined system. 
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Table 1. Geometric and structural properties 

Pr operty Coin}). 1 Comp. 2 Whole 
Length 15 5 20 

Mass/Length 6.775 6.775 6.775 

Rigidity, El 36458.3 36458.3 36458.3 

No. of modes 20 20 20 


Control Design 

The control design for component no. 1 is briefly 
described. As shown in Fig. 4, there are two 
control inputs available for this component. There 
is a torque wheel at the third of the length (x=5m) 
and there is a thruster at the middle of the beam 
(x=7.5m) to provide attitude and position control. 
Collocated with these actuators are attitude and 
position sensors, providing measurement for the 
feedback loop. Here, sensors and actuators dynamics 
are not considered. The control design was based on 
the Hm synthesis. The performance requirements 
were defined in terms of the weighted output sen- 
sitivity transfer function from disturbances to the 
measurement sensors. The weighting function was 
chosen as W p = ( ,"Sn^ To) for each channel to 
provide good disturbance rejection at low' to mid 
frequencies as well as to provide integral action to 
minimize steady-state errors. Uncertainty in the 
system model, in the form of unstructured input 
multiplicative, w r as used. The weighting function, 
used to scale or distribute the uncertainty, was 
chosen as W u = to emphasize model un- 

certainty in the mid-high frequency range. A 10- 
order model of the beam, which included the first 
three flexible modes, were used in the control design. 
The controller was synthesized using the ’’hinfsyn” 
routine of //-Analysis and Synthesis Toolbox. 8 The 
design converged at a 7 value of 0.942 (i.e., all 
the design requirements were met), and resulted in 
a 16th-order controller. Figure 5 provides a bode 
plot of the controller. This controller w T as later used 
in robustness analysis of the system in the docked 
configuration. 

Robustness Analysis 

Here, the two components w r ere assumed to be 
in a docked configuration. Each of the components 
was modeled using the first 20 modes (including two 
rigid-body modes and 18 flexible modes), resulting 
in a 40th-order state space model for each of the 
substructures. The 20 modes used w^ere deemed 
sufficient to provide a basis for describing the first 




Fig. 5 Bode Plot for the i /00 Controller for 
Component No. 1 

20 modes of the combined structure. The com- 
bined structure was also modeled globally (without 
synthesis) as a 20-m long beam using its first 20 
modes, resulting in a 40th-order state space model. 
This model was used for validation of the synthesis 
approach and is referred to as the truth model. 


Table 2. Natural frequencies 


Mode No. 

True 

Synthesized 

1 

0 

0 

2 

0 

0 

3 

0.0410 

0.0414 

4 

0.1131 

0.1158 

5 

0.2217 

0.2264 

6 

0.3665 

0.3676 

7 

0.5475 

0.5500 

8 

0.7647 

0.7809 

9 

1.0181 

1.0389 

10 

1.3077 

1.3108 

11 

1.6335 

1.6437 

12 

1.9955 

2.0425 

13 

2.3938 

2.4466 

14 

2.8282 

2.8344 

15 

3.2988 

3.3251 

16 

3.8056 

3.9045 

17 

4.3488 

4.4533 

18 

4.9278 

4.9389 

19 

5.5432 

5.6002 

20 

6.1948 

6.3764 


Using the procedure outlined in Eqs. (11)-(13), 
the dynamics of the two components were synthe- 
sized (combined). Table 2 compares the natural 
frequencies of the synthesized model vs. that of 
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sized System 

the validation (truth) model. Comparison of the 
frequencies indicate a good match between the two 
models. 

Now, both the synthesized model as well as the 
truth model were assessed against uncertainties in 
the frequencies and damping of the two components. 
Specifically, it was assumed that there is a five 
percent uncertainty associated with the frequency of 
the first mode in each component. Furthermore, it 
was assumed that there is a 15 percent uncertainty 
associated with the damping of the first mode in 
each component. The inherent structural damping 
ratio used for each of the component were chosen at 
0.5 percent. The damping ratio used for the truth 
model was based on the values obtained from the 
synthesized model. 

First, robust stability of the synthesized model 
was assessed. Following Eq. (23), the coeffi- 
cients Qi for the four uncertainty terms in the 
state matrices of the components, corresponding to 
frequency and damping in their first mode, were 
defined and factored according to Eq. (24). Note 
that the uncertainties were assumed to be distinct, 
real scalars. Using the system defined in Eq. 
(21) (with appropriate input/output matrices for 
the uncertainties) together with the controller 
designed for component no. 1, a robust stability 
analysis was performed with aid of the ”mu” routine 
of //-Analysis and Synthesis Toolbox. 8 A Mu plot 
of the closed-loop system is provided in Fig. 6, 
where it is observed that the structured singular 
value fi is less than one for all frequencies. Therefore, 
the synthesized system with the H^, controller is 
robustly stable. 

Next, robust stability of the truth model against 
the same uncertainties in the component models 
was assessed. The first task was to approximately 
translate the uncertainties from the component level 
to the system level. This was performed by trans- 
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Fig. 7 Robust Stability fi Plot for the Truth 
System 

forming the state matrix of the synthesized system 
to real modal form, with eigenvalues in increasing 
order on the diagonal. Then, Eq. (31) was used to 
determine the influence matrices at the system level, 
and the truth model given by Eq. (26) was used in 
robust stability analysis. Again, the H 0 0 controller 
designed for component no. 1 was the controller 
used in the robustness analysis of the closed-loop 
truth model. A Mu plot of the closed-loop system 
is provided in Fig. 7, where it is observed that the 
structured singular value // is less than one for all 
frequencies. Comparing Figs. 6 and 7, it is observed 
the Mu plots for the synthesized and truth systems 
match each other reasonably well. This validates 
the approximate characterization of the uncertainty 
at the system level from component level definitions. 

Concluding Remarks 

Robustness of flexible systems in the presence 
of model uncertainties at the component level was 
considered. Specifically, an approach for formulating 
robustness of flexible systems in the presence of fre- 
quency and damping uncertainties at the component 
level was presented. A controls-based algorithm for 
component mode synthesis was followed and modi- 
fied to avoid to superfluous modes that may come 
about in the dynamics. First, robustness of synthe- 
sized systems was formulated in a Linear Fractional 
Transformation framework. The formulation was 
then extended to deal with global (non-synthesized ) 
dynamic models with component-level uncertainties 
by projecting uncertainties from component levels to 
system level. A numerical example involving a two- 
dimensional simulated docking problem was worked 
out to demonstrate the feasibility of the proposed 
approach. The results demonstrated the effective- 
ness of the robustness formulation in transformed co- 
ordinates and the efficacy of the projection approach 
to translate uncertainties from their component-level 
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structures to system-level structures. 
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